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Abstract. - We present a computer simulation model that reproduces, event-by-event, the wave 
mechanical results of double-slit and two-beam interference experiments. The same model also 
simulates a one-to-one copy of a single-photon interference experiment with a Fresnel biprism 
(Jacques V. et ah, Eur. Phys. J. D, 35 (2005) 561). The model satisfies Einstein's criterion of 
local causality and is solely based on experimental facts, comprising the apparatuses used in the 
experiment and the observation of individual detector clicks. Our results prove that it is possible 
to give a particle-only description of single-photon double-slit experiments. 
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• Introduction. — In 1802, Young performed a double- 
slit experiment with light in order to resolve the ques- 
tion whether light was composed of particles, confirming 
Newton's particle picture of hght, or rather consisted of 
waves [1]. His experiment showed that the light emerg- 
ing from the slits produces a fringe pattern on the screen 
that is characteristic for interference, discrediting New- 
ton's corpuscular theory of light [1]. It took about hundred 
years until Einstein with his explanation of the photoelec- 
tric effect in terms of photons, somehow revived Newton's 
particle picture of light [2]. In 1924, de Broglie introduced 
the idea that also matter, not just light, can exhibit wave- 
like properties [3] . This idea has been confirmed in various 
double-slit experiments with massive objects such as elec- 
trons [4-7], neutrons [8,9], atoms [10,11] and molecules 
such as Cgo and C70 [12,13], all showing interference. The 
observation that matter and light exhibit both wave and 
particle character, depending on the circumstances under 
which the experiment is carried out, is reconciled by in- 
troducing the concept of particle- wave duality [14]. 

In most double-slit experiments, the interference pat- 
tern is built up by recording individual clicks of the detec- 
tors. In some of these experiments [5,6,15] it can be argued 
that at any time, there is only one object that travels from 
the source to the detector. Under these circumstances, the 
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real challenge is to explain how the detection of individual 
objects that do not interact with each other can give rise 
to the interference patterns that are being observed. 

According to Feynman, this phenomenon is "impossible, 
absolutely impossible to explain in any classical way and 
has in it the heart of quantum mechanics" [16]. Later, 
Feynman used the double-slit experiment as an example 
to argue that "far more fundamental was the discovery 
that in nature the laws of combining probabilities were not 
those of the classical probability theory of Laplace" [17]. It 
is known that the latter statement is incorrect as it results 
from an erroneous application of probability theory [18, 
19]. 

In this letter, we show that also Feynman's former state- 
ment needs to be revised: We present a simple computer 
algorithm that reproduces event-by-event, events being 
defined as clicks of a detector, just as in real double-slit 
experiments, the interference patterns that are usually as- 
sociated with wave behavior. We also demonstrate that 
our event-by-event simulation model reproduces the wave 
mechanical results of a recent single-photon interference 
experiment that employs a Fresnel biprism [15]. 

In our simulation model every essential component 
of the laboratory experiment such as the single-photon 
source, the slit, the Frensel biprism, and detector array 
has a counterpart in the algorithm. The data is analyzed 
by counting detection events, just as in Ref. [15]. The sim- 
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ulation model is solely based on experimental facts and 
satisfies Einstein's criterion of local causality. 

In a pictorial description of our simulation model, we 
may speak about "photons" generating the detection 
events. However, these so-called photons, as we will call 
them in the sequel, are elements of a model or theory for 
the real laboratory experiment only. The only experimen- 
tal facts are the settings of the various apparatuses and 
the detection events. What happens in between activat- 
ing the source and the registration of the detection events 
belongs to the domain of imagination. 

In the simulation model, the photons have which-way 
information, never have direct communication with each 
other and arrive one by one at a detector. Although the 
photons know exactly which route they followed, they nev- 
ertheless build up an interference pattern at the detec- 
tor, thereby contradicting the first part of Feynman's first 
statement, reproduced in the introduction. With our sim- 
ulation model, we demonstrate that it is possible to give 
a complete description of the double-slit experiment in 
terms of particles only. 

The event-based simulation approach that we describe 
in this letter is unconventional in that it does not re- 
quire knowledge of the wave amplitudes as obtained by 
solving the wave mechanical problem. Instead, the in- 
terference patterns are obtained through a simulation of 
locally causal, classical dynamical systems. Our approach 
provides a common-sense description of the experimental 
facts without invoking concepts from quantum theory such 
as the particle- wave duality [14]. 

In our simulation approach, we adopt the point of view 
that quantum theory has nothing to say about individual 
events [14]. Therefore, the fact that there exist event-by- 
event simulation algorithms that reproduce the results of 
quantum theory has no direct implications to the founda- 
tions of quantum theory: These algorithms describe the 
process of generating events at a level of detail that is 
outside the scope of what current quantum theory can de- 
scribe. The work presented here is not concerned with an 
interpretation or an extension of quantum theory. 

For phenomena that cannot (yet) be described by a de- 
ductive theory, it is common practice to use probabilis- 
tic models. Although Kolmogorov's probability theory 
provides a rigorous framework to formulate such mod- 
els, there are ample examples that illustrate how easy it 
is to make plausible assumptions that create all kinds of 
paradoxes, also for every-day problems that have no bear- 
ing on quantum theory at all [20-23]. Subtle mistakes 
such as dropping (some of the essential) conditions, like in 
the discussion of the double-slit experiment [18,19], mix- 
ing up the meaning of physical and statistical indepen- 
dence, changing one probability space for another dur- 
ing the cause of an argument can give rise to all kinds 
of paradoxes [18,23-26]. It seems that the interference 
phenomena that we focus on in this letter cannot be ex- 
plained within the framework of Kolmogorov's probability 
theory [27] . Our simulation approach does not rely on con- 
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Fig. 1: Schematic diagram of an interference experiment with a 
Fresnel biprism (FBP) [43]. S, Si, S2 denote the point source 
and its two virtual images, respectively. The grey area is the 
region in which an interference pattern can be observed. 



cepts of probability theory: We strictly stay in the domain 
of finite-digit arithmetic. 

Event-by-event simulation. — In our simulation 
approach, the photons are regarded as messengers that 
travel from a source to a detector. The messenger carries 
a message that may change as the messenger encounters 
another object on its path such as the Fresnel biprism 
for instance. The algorithms for updating the messages 
are designed in such a way that the collection of many 
messages, that is the collection of many detection events, 
reproduces the results of Maxwell's theory. The key point 
of these algorithms is that they define classical, dynamical 
systems that are adaptive. 

This letter builds on earlier work [28-37] that demon- 
strates that it is possible to simulate quantum phenomena 
on the level of individual events without invoking concepts 
of quantum theory or probability theory. Specifically, we 
have demonstrated that locally-connected networks of pro- 
cessing units with a primitive learning capability can sim- 
ulate event by event, the single-photon beam splitter and 
Mach-Zchnder interferometer experiments of Grangier et 
al. [38], Einstein-Podolsky- Rosen experiments with pho- 
tons [39-41], universal quantum computation [30,31], and 
Wheeler's delayed choice experiment of Jacques et al. [42]. 

In our earlier work, there was no need to simulate the 
detection process itself but, as we argue later, to simulate 
two-beam interference event by event, it is logically im- 
possible to reproduce the results of wave theory without 
introducing a model for the detectors. We show that the 
simplest algorithm that accounts for the essential features 
of a single-event detector allows us to reconstruct, event 
by event, the interference patterns that are described by 
wave theory. Incorporating this detector model in the sim- 
ulation models that we reported about in our earlier work 
does not change the conclusions of Refs. [28-37]. In this 
sense, the detector model adds a new, fully compatible, 
component to our collection of event-based algorithms. 

Wave theory. — Figure 1 shows a schematic dia- 
gram of a two-beam interference experiment with a Fresnel 
biprism [43]. A pencil of light, emitted by the source S, 
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Fig. 2: Schematic diagram of a double-slit experiment with two 
sources Si and S2 of width a, separated by a center-to-center 
distance d. The sources emit light according to a uniform in- 
tensity distribution and with a uniform angular distribution, 
P denoting the angle. The light is recorded by detectors D 
positioned on a semi-circle with radius X. 
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Fig. 3: Schematic diagram of a two-beam interference experi- 
ment with two line sources S\ and S2 having a spatial Gaus- 
sian profile. The sources are separated by a center-to-center 
distance d and emit light according to a uniform angular distri- 
bution, P denoting the angle. The light is detected by detectors 
D positioned at {X,y). 



is divided by refraction into two pencils [43] . Interference 
can be obtained in the region where both pencils overlap, 
denoted by the grey area in Fig. 1. As a Fresnel biprism 
consists of two equal prisms with small refraction angle 
and as the angular aperture of the pencils is small, we 
may neglect aberrations [43]. The system consisting of 
the source S and the Fresnel biprism can then be replaced 
by a system with two virtual sources Si and S2 [43], see 
Fig. 1. For such a system, it is straightforward to compute 
the intensity of the wave at the detector screen. 
We consider two cases: 

• The sources 6*1 and S2 are slits of width a, separated 
by a center-to-center distance d, see Fig. 2. Then, in 
the Fraunhofer regime, we have [43] 
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where A is a constant, k is the wave number, and 9 
denotes the angular position of the detector on the 
circular screen. 

The sources Si and 52 form a line source with a cur- 
rent distribution given by 



J{x,y) = 6{x) ^ 



-(a-sd/2)V2(T2 



(2) 



where a is the variance and d denotes the distance 
between the centre of the two sources. The intensity 
of the overlapping pencils reads 

^ , , ^ / , hyd (1 - h)kyd\ -b(y^+d^/4) 
I{y)^B{ cosh + cos ^ j e ^ 

(3) 

where S is a constant, h = k^a'^/{X^ + fc^cr^), and 
{X, y) are the coordinates of the detector (see Fig. 3). 
Closed-form expression Eq. (3) was obtained by as- 
suming that d <^ X and a ^ X . 



From Eqs. (1) and (3), it directly follows that the intensity 
distribution on the detector screen displays fringes that are 
characteristic for interference. 

Results of a time-resolved experiment that is a labora- 
tory realization of the interference experiment schemati- 
cally depicted in Fig. 3 are presented in Ref. [44]. 

Simulation model. — Looking at Figs. 2 and 3, com- 
mon sense leads to the conclusion that once we exclude 
the possibility that there is direct communication between 
photons, the fact that we observe individual events that 
form an interference pattern can only be due to the pres- 
ence and internal operation of the detector (we ignore the 
possibility that there is no common-sense explanation for 
this phenomenon). Therefore, it is worthwhile to consider 
the detection process in more detail. 

In its simplest form, a light detector consists of a mate- 
rial that can be ionized by light. The electric charges that 
result from the ionization process are then amplified and 
subsequently detected by appropriate electronic circuits. 
As usual, the interaction between the incident electromag- 
netic field E and the material takes the form P • E, where 
P is the polarization vector of the material [43]. Treat- 
ing this interaction in first-order perturbation theory, the 
detection probability is P{t) = J* ((E^(i') • K(i' - t") ■ 
I]{t")))dt'dt" where K{t' - t") is a memory kernel that 
contains information about the material only and ((.)) de- 
notes the average with respect to the initial state of the 
electromagnetic field [23]. Very sensitive photon detec- 
tors such as photomultipliers and avalanche diodes have 
an additional feature: They are trigger devices meaning 
that the generated signal depends on a threshold. 

From these general considerations, it is clear that a min- 
imal model should be able to account for the memory and 
the threshold behavior of real photon detectors. As it 
is our intention to perform event-based simulations, the 
model for the detector should, in addition to the two fea- 
tures mentioned earlier, operate on the basis of individual 
events. 



F. Jin et al. 



Messenger. In our simulation approach, we view each 
photon as a messenger. Each messenger carries a mes- 
sage, representing its time of flight. Compatibility with 
the macroscopic description (Maxwell's theory) demands 
that the encoding of the time of flight is modulo a dis- 
tance which, in Maxwell's theory, is the wavelength of the 
light. Thus, the message is conveniently encoded as a two- 
dimensional unit vector e,; = (eo,i,ei,j;) = (cos^i, sin0i), 
where (pi is the event-based equivalent of the phase of the 
electromagnetic wave and the subscript i > labels the 
messages. 

Source. A single-photon source is trivially realized in 
a simulation model in which the photons are viewed as 
messengers. We simply generate a message, wait until 
this message has been processed by the detector, and then 
generate the next message. This ensures that there can be 
no direct communication between the messages, implying 
that our simulation model (trivially) satisfies Einstein's 
criterion of local causality. 

Detector. A simple model for the detector that ac- 
counts for the three features mentioned earlier, contains 
an internal vector = {po,iiPi,i) with Euclidean norm 
less or equal than one. This vector is updated according 
to the rule 

P» = IVt-i + (1 - 7)ej, (4) 

where < 7 < 1. The machine generates a binary output 
signal Si using the threshold function 

S, = Q{pl,+vl^~n), (5) 
where < < 1. The total detector count is defined as 

M 

N = Y.S,, (6) 

i=l 

where AI is the total number of messages received. Thus, 
N counts the number of one's generated by the machine. 
Although not essential for our simulations, it is conve- 
nient to use pseudo-random numbers for r^. This will 
mimic the unpredictability of the detector signal. The 
parameter 7 controls the precision with which the mes- 
sage processing machine defined by Eq. (4) can represent 
a sequence of messages with the same and also controls 
the pace at which new messages affect the internal state 
of the machine [28]. The internal vector and 7 play 
the roles of the polarization vector P(t) and the mem- 
ory kernel K(t' — t"), respectively. Notice that the formal 
solution of Eq. (4) has the same mathematical structure 
as the constitutive equation P(i) = x(M)E(i — u)du 
in Maxwell's theory [43]. The threshold function of the 
real detector is implemented through Eq. (5). A detec- 
tor screen is just a collection of identical detectors and is 
modeled as such. Each detector has a predefined spatial 
window within which it accepts messages. 

It is not easy to study the behavior of the classical, 
dynamical system defined by Eqs. (4)-(6) by analytical 
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Fig. 4; Detector counts as a function of the angular (spatial) de- 
tector position (y) as obtained by event-by-event simulations 
of the interference experiment shown in Fig. 2 (Fig. 3). The 
open circles denote the event-based simulation results. The 
dashed lines are the results of wave theory (see Eqs. (1) and 
(3)). (a) The sources are slits of width a = A (A = 670 nm 
in all our simulations), separated by a distance d — 5A, see 
Fig. 2; (b) The sources Si and S2, separated by a distance 
d = 8A, emit single photons according to a Gaussian distribu- 
tion with variance a — X and mean d/2 and —d/2, respectively 
(see Fig. 3). 



methods but it is close to trivial to simulate the model on 
a computer. 

Note that the machine defined by Eq. (4) has barely 
enough memory to store the equivalent of one message. 
Thus, the machine derives its power from the way it pro- 
cesses successive messages, not from storing a lot of data. 
In particular, the machine does not know about M, a 
quantity that is unknown in real experiments. 

It is self-evident that the detector model defined by 
Eqs. (4) and (5) should not be regarded as a realistic model 
for say, a photomultiplicr. Our aim is to show that, in the 
spirit of Occam's razor, this is probably the simplest event- 
based model that can reproduce the interference patterns 
that we usually describe by wave theory. 

Simulation results. — First, we show that our event- 
by-event simulation model reproduces the wave mechan- 
ical results Eq. (1) of the double-slit experiment. Sec- 
ond, we simulate a two-beam interference experiment 
and demonstrate that the simulation data agrees with 
Eq. (3). Finally, we present the results for the simula- 
tion of the single-photon interference experiment with a 
Fresnel biprism [15], as depicted in Fig. 1. 

Double-slit experiment. As a first example, we con- 
sider sources that are slits of width a = A (A = 670 nm 
in all our simulations), separated by a distance d = 5 A, 
see Fig. 2. In Fig. 4(a), we present the simulation re- 
sults for a source-detector distance X = 0.05 mm and for 
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Fig. 5: Schematic diagram of the simulation setup of a single- 
photon experiment with a Fresnel biprism. The apex of the 
Fresnel biprism with summit angle a is positioned at (X',0). 
A line source with a current distribution given by Eq. (2) emits 
single photons from points S and with angles /3 chosen ran- 
domly from the interval [—a/2, a/2]. The detectors D posi- 
tioned at (X, y) count the photons. 



7 = 0.999. When a messenger (photon) travels from the 
source at (0,?/) to the detector screen positioned at X, it 
updates its own time of flight, or equivalently its phase 0^. 
This time of flight is calculated according to geometrical 
optics [43]. As the messenger hits a detector, the detector 
updates its internal state and decides whether to output 
a zero or a one. This process is repeated many times. 
The results of wave theory, as given by Eq. (1), are repre- 
sented by the dashed lines. Looking at Fig. 4(a), it is clear 
that there is excellent agreement between the event-based 
simulation and wave theory. 

Two-beam interference experiment. As a second ex- 
ample, we assume that the messengers leave either source 
5*1 or 5*2 from a position y that is distributed according to 
a Gaussian distribution with variance a and mean +d/2 or 
— d/2, respectively (see Fig. 3). The simulation results for 
a source-detector distance X = Q.l mm and for 7 0.999 
arc shown in Fig. 4(b). The dashed line is the correspond- 
ing result of wave theory, see Eq. (3). Also in this case, the 
agreement between wave theory and the event-by-event 
simulation is extremely good. 

Experiment with a Fresnel biprism. For simplicity, we 
assume that the source S is located in the Fresnel biprism. 
Then, the results do not depend on the dimension of the 
Fresnel biprism. Figure 5 shows the schematic represen- 
tation of the single-photon interference experiment that 
we simulate. Simulations with a Fresnel biprism of finite 
size yield results that differ qualitatively only (results not 
shown) . The time of flight of the ith message is calculated 
according to the rules of geometrical optics [43] . 

In the simulation, the angle of incidence f3 of the photons 
is selected randomly from the interval [—a/2, a/2], where 
a denotes the summit angle of the Fresnel biprism. A col- 
lection of representative simulation results for 7 = 0.999 
is shown in Fig. 6, together with the results as obtained 
from wave theory. Again, we find that there is excellent 
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Fig. 6: Detector counts as a function of the detector position y 
of the detector array positioned at X (see Fig. 5). The Fresnel 
biprism has an index of refraction n = 1.5631 and a summit an- 
gle a = 1°. Its apex is positioned at {X' , 0) with X' — 45 mm. 
The source emits single photons according to a Gaussian inten- 
sity distribution with variance a = 0.531 mm and wavelength 
A = 670 nm [15] . The open circles denote the event-based sim- 
ulation results. The dashed lines denote the results as obtained 
from wave theory. {&) X — X' = 7 mm; (b) X — X' = 15 mm; 
(c) X — X' — 55 mm. Thousand detectors where used to record 
the individual events. 



quantitative agreement between the event-by-event simu- 
lation data and wave theory. Furthermore, the simulation 
data presented in Fig. 6 is qualitatively very similar to the 
results reported in Ref. [15] (compare with Fig. 4(d) and 
Fig. 5(a)(b) of Ref. [15]). 

Conclusion. — In this letter, we have demonstrated 
that it is possible to give a particle-only description for 
single-photon interference experiments with a double-slit, 
two beams, and with a Fresnel biprism. Our event-by- 
event simulation model 

• reproduces the results from wave theory, 

• satisfies Einstein's criterion of local causality, 

• provides a pictorial description that is not in conflict 
with common sense. 

We do not exclude that there are other event-by-event 
algorithms that reproduce the interference patterns of 
wave theory. For instance, in the case of the single-electron 
experiment with the biprism [45], it may suffice to have 
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an adaptive machine handle the electron-biprim interac- 
tion without having adaptive machines modehng the de- 
tectors. We leave this topic for future research. 

We hope that our simulation results stimulate the design 
of new time-resolved single-photon experiments to test our 
particle-only model for interference. However, in order to 
eventually falsify our model, it would not suffice to simply 
publish raw experimental data and show that they do not 
agree in all details with the model described in this let- 
ter. Indeed, the model that we have presented is far from 
unique and it requires little imagination to see that one 
can construct similar but more sophisticated event-based 
models that can explain interference without making re- 
course to wave theory. 
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